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oscillator strengths of transitions from such states to the J = ground state. Hole energies of \aw-Z 
neonlike ions are evaluated. 
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I. INTRODUCTION 

A combined configuration-interaction (CI) many- 
body-perturbation-theory (MBPT) method, applied pre- 
viously to divalent atoms 0, is extended to particle- 
hole states of closed-shell atoms. After derivation of 
CI+MBPT expressions for particle-hole states, we will 
apply the theory to calculations of energies and electric- 
dipole transition probabilities for neon. 

For neon, many accurate measurements of transition 
rates are available, providing important tests of theory. 
Reciprocally, the theory might help resolve existing dis- 
crepancies among oscillator strengths (/-values) for tran- 
sitions from the ground state to several excited states, for 
which experiments disagree. There is also a certain de- 
ficiency in existing ab-initio theories in neon, for which 
discrepancies among many measurements and theoret- 
ical calculations are unsettled. For example, the only 
other elaborate ab-initio calculations (Avgoustoglou and 

Beck |]) give an oscillator strength for the 1 2^3^3,51/2 1 

neon state larger than most experimental values by more 
than two standard deviations. Extensive calculations 
performed by Hibbert et al. || for many transition 
rates along the neon isoelectronic sequence use a general 
configuration- interaction code (CIV3) 0. The calcula- 
tions utilize parametric adjustments with measured fine 
structures, but do not completely agree with experiments 
in neon and have an accuracy similar to other semicmpir- 
ical calculations of Seaton M. However, the two calcula- 
tions disagree with each other for several transitions. We 
hope that our calculations may help to understand better 
the theoretical problems in neon and provide guidance for 
the analysis of experimental data. 

Some possible applications of the present CI+MBPT 
method include the study of neonlike ions, Ne I - Si V, S 
VII, Ar IX, Ca XI, and Fe XVII that have astrophysical 



interest and have been included in the Opacity Project 
(Seaton Q ) . The transition data in neon and other noble 
gases are also used in plasma physics, and in studying dis- 
charges that find many industrial applications in lamps 
and gas lasers. The methods presented here might be also 
used for improving the accuracy of MBPT or for extend- 
ing CI+MBPT to more complicated open-shell atoms. 

The principal theoretical difficulty arises from the sen- 
sitivity of transition amplitudes to the interaction be- 
tween closely spaced fine-structure components. Al- 
though it is possible to obtain energies which are rea- 
sonably precise on an absolute scale using coupled- 
cluster methods (Ilyabaev and Kaldor Q| ) , accurate fine- 
structure splittings seem very difficult to obtain without 
semiempirical adjustments. This is why semiempirical 
approaches, which have fine-structure intervals carefully 
adjusted, are more successful in neon than are ab-initio 
calculations. However, as we will demonstrate in this 
paper, CI calculations corrected with MBPT are also 
capable of accurately predicting fine-structure splittings 
and, consequently, transition amplitudes. In this paper, 
we will demonstrate the excellent precision of CI plus 
second-order MBPT. Third-order corrections, for which 
numerical codes already exist ||, can also be included, 
providing even further improvement in accuracy. 

In the following section, we use the effective Hamil- 
tonian formalism and particle-hole single-double cou- 
pled equations to derive expressions for the second-order 
Hamiltonian matrix of the CI+MBPT method. In the fi- 
nal expressions, we present a quite accurate new MBPT 
that can predict energies of hole states and can describe 
appropriately the interactions in particle- hole atoms. 
The accuracy of hole energies obtained with the new 
MBPT will be illustrated for neon and low-i? neon-like 
ions. Our CI+MBPT energies and /-values for many 
states of neon are tabulated. Their agreement with ex- 
periment and other theories are shown. 
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II. CI+MBPT METHOD 

The accuracy of the Raylcigh-Schrodinger variant of 
second-order MBPT given in |9[] is insufficient for pur- 
pose, so that more accurate single-double equations must 
be used. The formulas for the correlation operator and 
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a system of coupled equations for the correlation coef- 
ficients are given in we follow the notation of 
in the the paragraphs below. Under certain conditions, 
those equations can be further simplified and rewritten 
in the following form: 
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In the second equation of this set, the term cjbcbcXbc * s 
subtracted from both sides of this equation to make the 
right-hand side small. Since large random-phase approxi- 
mation (RPA) corrections in the particle-hole CI+MBPT 
are treated by CI, the quantities W" a ? entering this set 
of equations on the right-hand side in Ref. JhJ are small 
and have been neglected here. The concern might be 
raised for the correlation coefficients Xvl anc ^ Xvi which 
generally would have small factors (e v + Eb — £ a — £s) or 
(e v — ey) in front. However, for the large CI model space, 
energies of the core- virtual orbitals bs are well separated 
from the energies of the valence-hole orbitals av. The 
quantities R in zero approximation can be set to: 
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to obtain the first-order effective Hamiltonian, 

Kw,va = (St, " £a)6 v > v Sa'a + (3) 

and the correlation coefficients x- Here we define the 
first-order correction Hy aa , v = A, v > v 5 a > a + g v > a a>v to the 
effective Hamiltonian. For faster convergence of CI and 
for subtraction of the dominant monopole contributions 
in RPA diagrams, a V( N -^ Hartree-Fock (HF) model 
potential for which A nm = g nam a, A na = A an = A ab = 
is introduced. 

Further improvement of accuracy can be achieved 
through iterations. After one iteration we obtain the 
second-order contribution to the effective Hamiltonian, 

H^Lv - 5Ri8 a>a + 5Ri5v>v + SKI', (4) 



Note that in the last equation we have extended the 
single-double method. The last term entering 8R™ in 
the single-double formalism would normally not contain 
Hadad m the denominator. However, if we do not modify 
this denominator, we find that in the third-order MBPT, 
large terms proportional to g a dad will appear leading to 
a decrease in accuracy. A physical reason for modify- 
ing the denominator of this term is that the process de- 
scribed by this term contains two holes in the intermedi- 
ate states with large interaction energy. This interaction 
should be treated nonperturbatively, for example, by in- 
clusion of g a dad into the denominator as we have done on 
the basis of the single-double equations in other terms. 
Finally, this term is almost equal to the seventh term 
(they are complex conjugates and their Goldstone dia- 
grams are related by a reflection through a horizontal 
axis), and for convenience they are set equal in numer- 
ical calculations. The angular reduction for 5R%% can 
be easily obtained using the second-order particle-hole 
formulas given in Ref. ||. 
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III. A SOLUTION OF THE HOLE-ENERGY 
PROBLEM 

A. Breit corrections 

Apart from Coulomb correlation corrections, the Breit 
magnetic interaction is also important in neon and the 
isoelectronic ions. The breakdown of various Coulomb 
and relativistic contributions to the energy of 3s states 
of neon are given in Ref. |l(|. Breit corrections cancel, 
but for higher excited states they may not. Hence, to 
improve the accuracy of fine-structure splittings, we in- 
clude the Hartree-Fock hole Breit correction -B£i F ' in our 
calculations, 



TABLE I: A comparison of theoretical and experimental hole 
energies and the 2p3/2-2p!/ 2 fine-structure intervals for neon 
and neon-like ions. All energies are in cm 





Ne 


Na+ 


Mg + ' 2 


Al+ y 


Si+ 4 


2p 3 /2 Th. 
2p 3 /2 Exp. 
Difference 


172434 
173930 
1496 


380443 
381390 
947 


645951 
646402 
451 


967531 
967804 
273 


1344344 
1345070 
726 


2pi/ 2 Th. 
2pi/ 2 Exp. 
Difference 


173218 
174 710 
1492 


381816 
382756 
940 


648196 
648631 
435 


970997 
971246 
249 


1349449 
1350160 
711 


2p 3 /2-2pi/2, Th. 
2p 3 /2-2pi/2, Exp. 
Difference 


784 
780 
-4 


1373 
1366 
-7 


2245 
2229 
-16 


3466 
3442 
-24 


5090 
5105 
-15 



D aa " ac 



(8) 



We have checked that the first-order corrections to 
the energies of J = 2 and J — 1 states given in Table 1 
of Ref. [ fl0| agree with our fiia F) contributions, 0.00062 
and 0.00090 a.u., for 2p 3 / 2 and 2px/2 states, respectively. 
We omit the small frequency-dependent Breit, quantum- 
electrodynamic, reduced-mass, and mass-polarization 
corrections. Small as they are, those corrections are fur- 
ther reduced after subtraction for the fine-structure in- 
tervals. More careful treatment of relativistic corrections 
is needed in calculations of high-Z neon-like ions. 



B. Calculations of hole energies for neonlike ions 



TABLE II: A comparison with experiment of CI+MBPT ener- 
gies referenced to the ground state and given in atomic units. 
An almost constant shift is subtracted in the fifth column 
to demonstrate excellent agreement for relative positions of 
levels 



Level 

P3/2 3S 

P3~/2 4S 
Pl/2 4S 
P3/2 3d 

p~J 2 5s 



Experiment CI+MBPT A A - 0.0069 



0.6126 


0.6048 


0.0078 


0.0009 


0.6192 


0.6116 


0.0076 


0.0007 


0.7235 


0.7166 


0.0070 


0.0001 


0.7269 


0.7200 


0.0069 


0.0000 


0.7360 


0.7289 


0.0070 


0.0001 


0.7365 


0.7294 


0.0071 


0.0002 


0.7401 


0.7330 


0.0071 


0.0002 


0.7560 


0.7491 


0.0069 


0.0000 


0.7593 


0.7525 


0.0069 


0.0000 



Since we propose a new variant of the MBPT expan- 
sion, we would like first to demonstrate that this expan- 
sion is convergent for hole states. The theoretical hole 
energies shown in Table [Shave been obtained in the yW 
HF potential using Eq. (Q) for SR^ to calculate second- 
order corrections. The extra term in the denominator is 
important and is necessary for convergence of the per- 
turbation expansion. Experimental hole energies in the 
National Institute of Standards and Technology (NIST) 
database Ref. Jll| are found as the limit energies for the 
neon isoelectronic sequence. For neutral neon only one 
limit, the p 3 / 2 energy is given in NIST (H]]. The 2px/ 2 - 
2p3/2 splitting 780.4269(36) cm -1 has been measured in 
Ref. pj, and using this value we find the experimental 
Pi/2 energy. Table ffl demonstrates the good agreement of 
our theoretical P3/2, P1/2 energies as well as the same fine 
structure interval for neon- like ions. Our fine structure 
interval, whose correctness is crucial for transition ampli- 
tude calculations, differs from experiment just by about 
10 cm -1 . Note that the HF value 187175 cm -1 for the 
2p3/2 state is 8.5% higher than the experimental value 
173930 cm -1 , and, after adding correlation corrections, 
we obtain improvement by a factor of ten. For the fine 
structure, the HF value 1001cm -1 disagrees even more, 
by 28%. If we use Rayleigh-Schrodinger perturbation 



theory, the corrections are twice as large as our results, 
and the agreement with experiment does not improve. 



IV. NEON ENERGIES AND OSCILLATOR 
STRENGTHS OF J=l ODD STATES 

To test the accuracy of the CI+MBPT method, we 
first calculated energies of several lowest odd J=l neon 
states, Table [n| The number of configurations in CI was 
chosen to be 52. The order of eigenstates obtained in 
CI+MBPT is the same as the order of the experimental 
levels. We abbreviate long NIST designations since the 
levels are uniquely specified by energy or by order. 

The pure ab-initio energies differ from experimental 
energies by 0.0069 a.u., but after subtraction of the sys- 
tematic shift (which does not make much difference in 
transition calculations), the agreement is at the level of 
0.0001 a.u. for almost all states. Therefore, we consider 
the accuracy of CI+MBPT adequate for correct predic- 
tion of level mixing and oscillator strengths. For the 3s 
states, agreement with experiment for the fine structure 
interval is much better than that obtained by Avgous- 
toglou et al. 0.0002 versus 0.0012 a.u.; a possible 
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TABLE III: Our CI+MBPT oscillator strengths for the 
ground to excited state transitions in neon compared with 
average experimental values (3rd and 4th columns) and those 
obtained with the best semiempirical theories ^, q, |l| 



Levels CI+MBPT <r-avr mean 



p~ /2 3s 

P 3 ~/2 4s 
Pi/2 4s 
Pa/2 3d 

P3/\ 3d 

p 3/ , 2 5s 

Pl/2 5S 



0.0102 
0.1459 
0.0131 
0.0181 
0.0066 
0.0130 
0.0069 
0.0068 
0.0053 



0.0099 
0.1549 
0.0122 
0.0170 

0.0187 
0.0067 
0.0064 
0.0043 



0.0107 
0.1487 
0.123 
0.016 

0.0199 
0.0069 
0.0066 
0.0044 



Ref. |j Ref. || Rcf. 



0.0126 
0.1680 
0.0152 
0.0193 
0.0056 
0.0167 
0.0086 
0.0073 
0.0050 



0.0106 
0.1410 
0.0124 
0.0160 
0.0045 
0.0131 
0.0064 
0.0060 
0.0043 



0.0123 
0.1607 



0.0047 
0.0117 
0.0055 



TABLE IV: References for experimental data shown in Fig. |l| 



Obs. Reference 



Year 



/ 



1 Kuhn et al. |LJ 

2 Lawrence and Liszt jlf| 

3 Geiger @ 

4 Kernahan et al. 

5 Kazantsev and Chaika 116 

6 Knystautas and Drquin IE 91 

7 Bhaskar and Luriq_ 

8 Westerveld et al. 

9 Aleksandrov et al 

10 Chornay et al 

11 Tsurubuchi et al 

12 Chan et al. Q 

13 Ligtenberg et al 

14 Suzuki et al. ]27 

15 Curtis et al. [E8[ 

16 Gibson and Risley 

17 Zhong et al. [El 




1967 
1969 
1970 
1971 
1971 
1974 
1976 
1979 
1983 
1984 
1990 
1992 
1994 
1994 
1995 
1995 
1997 



0.01200 0. 
0.00780 0. 
0.00900 0. 
0.00840 0. 
0.01380 0. 
0.00780 0. 
0.01220 0. 
0.01090 0. 
0.01200 0. 
0.01200 
0.01220 
0.01180 
0.01070 
0.01060 0. 
0.00840 0. 
0.01095 0. 
0.01240 0. 



00200 
00040 
00200 
00070 
00080 
00080 
00090 
00080 
00300 
00400 
00060 
00060 
00030 
00140 
00030 
00032 
00380 



explanation for this could be that single-double equa- 
tions miss important corrections which we included by 
modifying the denominators. In Ref. pf]| , however, the 
systematic shift is small. 

Finally, we present our CI+MBPT oscillator strengths 
in neon. After diagonalization of the second-order effec- 
tive Hamiltonian, we obtain wave functions in the form 
of expansion coefficients in the CI space and use them 
to calculate oscillator strengths. Size-consistent formu- 
las for dipole matrix elements for transitions decaying 
into the ground state are provided in Ref. j| , where the 
absorption oscillator strength / is also defined. We give 
in this table ab-initio values of the oscillator strengths 
/. The dominant part of the RPA corrections is included 
at the level of CI. Small normalization corrections are 
omitted. 

Many experiments have disagreements in oscillator 
strengths far exceeding the cited errors (see Fig. |l| and 
Table IV): hence, for comparison, we give in Table § 
two statistical averages: the first is a weighted accord- 



0.017 
0.016 
0.015 
=! 0.014 

CO 

of 0.013 
si 

g> 0.012 

CD 

4= 0.011 

O) 

o 0.010 
f3 

= 0.009 
o 

O 0.008 
0.007 
0.006 




~i — 1 — i — 1 — i — 1 — i — 1 — r~ 

10 12 14 16 18 



2 4 6 

Observations in chronological order 



FIG. 1: Comparison with experiment for oscillator strengths 
of the [p^ 2 3s]i state of neon 



ing to cited standard deviations and the second is an 
unweighted average. For the 3s levels, the experimental 
data compiled in Ref. Q and for the higher excited levels 
in Ref. [ pOfl have been included in the averaging. Average 
values obtained here are not necessarily the most accu- 
rate, but they serve well for comparison and for a test of 
our probably less accurate calculated values. 

A more careful analysis of experimental techniques to 
exclude systematic errors, which are definitely present, is 
necessary; our values can provide some guidance. For 
p^ 2 3d states, since the energy separation of the two 
states is small, experiments give the sum of the two oscil- 
lator strengths, and the value 0.0196 rather than 0.0130 
should be compared with the experimental values 0.0187 
(0.0199). In this table, we also compare our theory with 
other semiempirical theories. Surprisingly, early calcu- 
lations by Aleksandrov et al. |p|] agree well with our 
calculations. A fair agreement, considering the high sen- 
sitivity of these transitions to correlation correction, is 
also obtained with the other theories in the table. 



V. CONCLUSIONS 

In this paper, we have introduced CI+MBPT theory 
for particle-hole states of closed-shell atoms. A difficulty 
that the hole energy has poor convergence is overcome 
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with modifications of denominators in MBPT. Good pre- 
cision for hole states and for particle-hole states is illus- 
trated for many energy levels of neon. Apart from en- 
ergies, our theory is tested in calculations of oscillator 
strengths. Agreement with averaged experimental values 
is achieved. 
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